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Abstract 

I review progress on selected issues connected with generalized parton distributions. Topics 
range from the description of hard exclusive reactions to the spatial distribution of quarks in the 
nucleon and the contribution of their orbital angular momentum to the nucleon spin. 


1 Introduction 

An outstanding task in QCD is to understand hadron structure at the level of quarks and gluons. A 
wealth of information about this structure is contained in generalized parton distributions (GPDs). 
They are related both to the conventional parton densities and to elastic form factors—quantities 
which have played a crucial role in the understanding of strong interactions, and which at first sight 
are of very different nature. Recent reviews on GPDs can be found in mm- 

GPDs are defined through matrix elements {p'\0\p) between hadron states 1^) and \p), with non¬ 
local operators O constructed from quark and gluon fields. As sketched in Fig. GPDs depend on 
several kinematical variables, namely on the longitudinal momentum fractions x + ^ and a: — ^ of the 
partons and on the invariant momentum transfer t = {p — p')"^ ■ For unpolarized quarks there are two 
distributions H‘>(x,^,t) and E'^(x,^,t). The former is diagonal in the proton helicity, whereas the 
latter describes proton helicity flip. For p = p' and equal proton helicities one recovers the diagonal 
matrix element parameterized by usual quark and antiquark densities, so that H‘^(x,0,0) = q{x) and 
H‘>{—x,0,0) = —q{x) for x > 0. Taking Mellin moments of in x, one obtains matrix elements of 
local operators, which are parameterized by form factors. In particular, the lowest moments give the 
well-known Dirac and Pauli form factors: 

'^eq f dxH'^{x,C,t) = Fi{t), '^eq f dxE'^{x,^,t) = F 2 {t), (1) 

q q •' 

where denotes the fractional quark charge. The ^ independence of the integrals is a consequence 
of Lorentz invariance. Of great interest is the next highest moment 

J dx x[R''^(x, t) + F^(x, t)] = 2 (2) 

because J'^(O) gives the total angular momentum carried by quarks and antiquarks of a given flavor, 
including both their helicity and their orbital angular momentum. Belonging to local operators, the 
Mellin moments of GPDs are well suited for evaluation in lattice QGD (see below). 
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Figure 1: a: Variables in a GPD. The momentum fractions x and ^ refer to the average hadron 
momentum \{p + p')- b: A graph for the amplitude of deeply virtual Compton scattering (DVCS). 


Factorization theorems state that GPDs appear in the scattering amplitudes of suitable hard 
exclusive processes. Figure^^ shows an example graph for deeply virtual Compton scattering, 7 *p ^ 
7 P, where a hard scale is provided by the initial photon virtuality. A large number of reaction 
channels can be accessed in hard exclusive meson production, such as 7 *p —> pp or 7 p —> 
where in addition to quark or gluon GPDs the quark-antiquark distribution amplitudes of the meson 
appear as non-perturbative input (see Fig. n. The longitudinal momentum transfer ^ is hxed by the 
process kinematics (for DVCS and light meson production one simply has ^ = xs/(2 —in terms of 
the usual Bjorken variable). In contrast, x is a loop variable, and ^ only gives the typical size of x in 
the convolution of the GPD with the hard-scattering amplitude. The hard-scattering subprocesses for 
Compton scattering, electroproduction of light mesons, and photoproduction of heavy quarkonium 
are fully calculated at next-to-leading order (NLO) in Ug, and partial results for DVCS at NNLO have 
just appeared [3]. Whereas ag corrections are generally found to be moderate for Compton scattering, 
they can be substantial for meson production |^, and more detailed studies will be needed to gain 
quantitative theoretical control over these channels. 

In addition to the variables already discussed, GPDs depend on the scale at which the partons 
are resolved, given by the hard scale in the physical process. The evolution equations for GPDs 
interpolate between those for ordinary parton densities and those for meson distribution amplitudes. 
The evolution kernels are known to NLO in ag, but it is only recently and only to LO accuracy that 
one has an explicit solution of the evolution equations [31 iOj- Generalizing the well-known Mellin 
transform technique employed for the ordinary parton densities, this solution should be useful both 
for efficient numerical implementation and for analytic considerations, for instance in the limit of 
small X. 

A typical strategy for modeling GPDs is to take an ordinary parton density as input and to generate 
a ^ dependence consistent with Lorentz invariance relations, which state that Mellin moments in x 



Figure 2: Example graphs for exclusive meson production. The large blobs represent GPDs and the 
small ones the distribution amplitudes of the produced meson. 
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of GPDs must be polynomials in ^ of a specified degree. Different ansatze are used to this end, the 
most common ones being based on so-called double distributions [7j or on the evolution of GPDs at 
small X and ^ jH]. Different strategies based on moments have recently been proposed in [Hj and in [Hj. 
With the Compton amplitude calculated at LO in Og, the model in provides a good description of 
the data dnum for DVCS at small xb- The t dependence of GPDs will be discussed below. 


2 Vector meson production 


An important feature of exclusive vector meson production is that quark and gluon distributions 
appear at the same order in ag, whereas in DVCS (just as in inclusive deep inelastic scattering) the 
gluon contribution is suppressed by ag relative to quarks. Schematically, one has 




+ u) + + d) + -^g 


(3) 


where u + u, d + d, s -1- s, and g represent the convolutions of the relevant GPDs with the hard- 
scattering kernels (which are identical for quarks and gluons at LO). Taking the ordinary parton 
densities as a rough guide for the relative size of terms in ©, one expects that, even in kinematics of 
the fixed-target experiments at HERMES and Jefferson Lab, gluons provide a substantial contribution 
to the production amplitude. This expectation is supported by preliminary HERMES data for the 
ratio of (j) and production cross sections m- The same conclusion has been reached in an explicit 
calculation with GPD models based on double distributions m- In the same analysis, substantial 
changes in the calculated cross sections at ~ 4 GeV^ and xs ~ 0.1 were found when taking different 
parameterizations of the gluon density as input to the GPD model. This reflects on one hand the 
present uncertainty in our knowledge even of the ordinary gluon distribution, and on the other hand 
the sensitivity of meson production to the generalized gluon distribution. Very similar findings were 
made in the study ^3] of exclusive production at collider energies: the high-precision data |15j 
clearly disfavor a number of parameterizations of the gluon density in conjunction with the Shuvaev 
model (HI for the ^ dependence of the gluon GPD. 

It has long been known that in exclusive meson production the leading-twist approximation, 
which gives the leading term in a 1/Q expansion, overshoots the data by factors of several for photon 
virtualities below 5 GeV^. Closer analysis reveals that important corrections to this approximation 
can be ascribed to the effect on the hard-scattering amplitude of the intrinsic transverse momentum 
of quarks in the produced meson. A recent analysis modeling this effect finds good agreement 
with p^ and (j) production data at high energies from HI and ZEUS. 


3 Impact parameter distributions 

An important feature of GPDs is that they contain information about the spatial distribution of 
quarks and gluons in the nucleon. This becomes explicit in the impact parameter representation im 
(for different approaches see CHI and d). To introduce this representation, let us form wave packets 

|P+.i>)=/(0e-'»'|p+,P> (4) 

from momentum eigenstates Ip"'', p), where the plus-momentum p+ = (p°-|-p^)/\/2 simply becomes the 
longitudinal momentum (up to a factor \/2) in a frame where the proton moves fast. The conjugate 
variable to the transverse momentum p = (p^,p^) is called impact parameter b and gives the position 
of the wave packet in the transverse plane. Indeed, |p“'",b) is an eigenstate of a suitably defined 
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transverse position operator, i.e., a relativistic particle can be localized exactly in two dimensions 
(without the ambiguities at the order of a Compton wavelength occurring when one attempts to 
localize a particle in all three dimensions). A more detailed analysis shows that b is the “center of 
momentum” of the partons in the proton, given by a weighted average b = pf / Si pf ™ terms 
of their plus-momenta and transverse positions. The center of momentum is related by Noether’s 
theorem to a class of Lorentz transformations called “transverse boosts”, in analogy to the relation 
between the center of mass and Galilean transformations in nonrelativistic mechanics. 

Forming matrix elements from impact parameter states and the operators defining generalized 
parton distributions in momentum space, one obtains Fourier transforms of these distributions. For 
^ = 0 one finds that 

q{x, b^) = J-^ 0, - A^) (5) 

is the density of quarks with longitudinal momentum fraction x and transverse distance b from 
the center of momentum of the proton. Integrating q{x,b‘^) over b one recovers the usual quark 
distribution. For ^ 7 ^ 0 one no longer has a probability interpretation because the two momentum 
fractions in Fig. ^ are different, but b still describes the distribution of the struck parton in the 
transverse plane. According to the discussion in the introduction, the combined ^ and t dependence 
of hard exclusive processes thus gives information about the impact parameter distribution of partons 
with longitudinal momentum fraction of order Precise measurements are in particular available for 
the t dependence of J/T production ^ 3 , giving access to the spatial distribution of small-x gluons. 
HI has published a first result for the t dependence of DVCS m, which at the small of the 
measurement is sensitive to a combination of sea quark and gluon distributions. 

Like ordinary parton densities, the impact parameter distributions q{x,b‘^) and g{x,b'^) for quarks 
and gluons depend on the resolution scale g. Their scale evolution is local in b and described by the 
usual DGLAP equations. As a consequence, the b dependence of the distributions at given x also 
changes with /i. A useful quantity to characterize the b distribution is the average squared impact 
parameter, defined as 


( 62 ) 


/ (fb 6^ q{x, b^) 
f (Pb q{x, 6^) 


4^ log H‘>{x,0,t) 


( 6 ) 


for quarks and in analogy for gluons. It is straightforward to obtain the evolution equations for {b‘^)x 
from those for the impact parameter distributions m- 

When a quark takes most of the proton momentum, its impact parameter will tend to coincide 
with the center of momentum of the proton as a whole. For x ^ 1 one therefore expects a narrow 
distribution in 6 , or equivalently a flat t dependence of GPDs in momentum space. An estimate for 
the overall transverse size of the proton in that limit is given by the transverse distance 6/(1 — x) 
between the struck quark and the center of momentum of the spectator partons, as shown in Fig. |S1 
It seems plausible to assume that this distance remains finite due to confinement EH, so that the 
average squared impact parameter of partons should vanish like {b‘^)x ~ (1 — x)^ for x —> 1. In the 
opposite limit of small x, the phenomenology of high-energy hadronic reactions suggests a behavior 
like of momentum-space GPDs PEO]. According to m this translates into a logarithmic 

growth of the average impact parameter as x becomes small, {b‘^)x ^ B + a'log(l/x). 

Indirect information on impact parameter distributions can be obtained by using the sum rules 
o and the extensive and precise data on the nucleon form factors EOlllll- For the Dirac form factor 
Fi{t) this requires an ansatz for the functional dependence of H'^{x,0,t), which can be restricted to 
its valence part, since the electromagnetic current is only sensitive to the difference of quark and 
antiquark distributions. A crucial result of such studies is the rapid decrease of the average impact 
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Figure 3: Left: Three-quark configuration with one fast quark in a proton. The thick line denotes the 
center of momentum of the proton and the dashed line the center of momentum of the two spectator 
quarks. Right: Average impact parameter of the valence u quark distribution (lower curve) 

and the associated distance (1 — between the struck parton and the center of momentum 

of the spectators (upper curve), as estimated in pU] . 


parameter with x over the entire x range, as illustrated for the distribution u{x, b‘^) — u{x, 6 ^) in Fig.|21 
In momentum space, this corresponds to a t dependence that becomes less steep with increasing x. 
This finding is confirmed by calculations in lattice QCD, where a clear decrease of the t slope is seen 
for moments / dx x'^H‘^(x,0,t) with increasing power n j23j . 

A strong correlation between the transverse distribution of partons and their momentum fraction 
is not only interesting from the perspective of hadron structure, but also has practical consequences 
for high-energy hadron-hadron collisions m Consider the production of a high-mass system (e.g. a 
dijet or a heavy particle). For the inclusive production cross section, the distribution of the colliding 
partons in impact parameter is not important: only the parton distributions integrated over impact 
parameters are relevant according to standard hard-scattering factorization (see Fig. |1^). There 
can however be additional interactions in the same collision, especially at the high energies of the 
Tevatron or the LHC, as shown in Fig. Their effect cancels in sufficiently inclusive observables, 
but it does affect the event characteristics and can hence be quite relevant in practice. In this case, 
the impact parameter distribution of partons does matter: the production of a heavy system requires 
large momentum fractions for the colliding partons. A narrow impact parameter distribution for these 
partons forces the collision to be more central, which in turn increases the probability for multiple 
parton collisions in the event. 



Figure 4: a: Graph with a single hard interaction in a hadron-hadron collision. The impact pa¬ 
rameters bi and 62 are integrated over independently. b: Graph with a primary and a secondary 
interaction. 
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Figure 5: A transition where the proton helicity is flipped but the quark helicity conserved. The he- 
licity mismatch is compensated by one unit of orbital angular momentum, so that angular momentum 
is conserved. 

4 Spin and orbital angular momentum 

The proton helicity-flip distributions have a density interpretation at ^ = 0, similar to the distri¬ 
butions H'J discussed so far. To see this one changes basis from proton helicity states |t), |i) to states 
|A±) = (ID ± II) )/'/2 polarized along the positive or negative x axis. In impact parameter space 
one then obtains the density 

uy n 

qx{x,b) = q{x,b^) -('^) 

m ob^ 

of unpolarized quarks in a proton polarized along the positive x direction, where the Fourier transform 

&') = / ^ E^x, 0, -A^) (8) 

is dehned in analogy to © and m denotes the nucleon mass. The impact parameter distribution 
of quarks in a transversely polarized proton is thus shifted in the direction perpendicular to the 
polarization. This shift must be substantial at least for some x and b: due to the sum rule o 
the average f dx f Sbe^{x,b‘^) = is given by the magnetic moments of proton and neutron as 
~ 1.67 and ~ —2.03 and hence quite large. A connection has been proposed in [2^1 between 
the anisotropy 0 in the spatial distribution of quarks and the Sivers effect, which is an anisotropic 
transverse momentum distribution of quarks in a transversely polarized proton. The Sivers effect has 
indeed been observed experimentally m- Similar asymmetries can be discussed for the distribution 
of transversely polarized quarks in an unpolarized proton and first results in lattice QCD indicate 
that the corresponding anisotropy of the impact parameter distribution is appreciable m- 

The interpretation of © as a density implies positivity conditions for e^{x,b‘^) [29| . In momentum 
space one has a simple bound 

\E'^{x, 0,0)1 < q{x) m^J{b^)^ , (9) 

with more stringent inequalities involving also quark helicity dependent distributions on the right- 
hand side. According to the behavior of {b‘^)x discussed in the previous section, this restricts E'^ quite 
severely at larger values of x. The distribution E'^ involves one nnit of orbital angular momentum 
since in the associated matrix elements the proton helicity is flipped but the quark helicity conserved 
(see Fig. El). The bound thus limits the amount of orbital angular momentum that can be carried 
by quarks m, and is especially strong at large x. 

Further information on the helicity-flip distribution can be obtained from the sum rule for the 
Pauli form factor in O- For this one needs an ansatz for the valence part of E'^{x,0,t), which can 
be made in analogy to the case of H'^(x,0,t) and the Dirac form factor [201122] . The corresponding 
fits are less well constrained, because in contrast to the forward limit E'^{x, 0,0) is not known and 
needs to be parameterized in addition to t dependence of E’^. It tnrns out that the positivity bound 
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© and its stronger versions provide valuable restrictions of the allowed parameter space. One can 
use these restrictions and estimate the total angular momentum carried by valence quarks, taking 
the analog of Ji’s sum rule © for the difference of quark and antiquark contributions. The result 
obtained in pni OT] is 2 = 0.39 0.46 and 2J^ = —0.04 0.04 at scale ^ = 2 GeV. Subtracting 

the quark helicity parts, one gets rather large numbers 2L“ = —(0.47 0.54) and 2L() = 0.30 0.38 

for the orbital angular momentum carried by u and d valence quarks. The corresponding estimate 
for the isovector combination 2(L“ — L^) = —(0.77 0.92) is large, whereas the isoscalar combination 

2(L“ + L(() = —(0.11 0.22) is uncertain and quite small. Lattice calculations by the QCDSF 

Collaboration obtain 2(L“ — L^) = 0.90 ± 0.12 and compatible with zero within errors [32j. 

Given the very different sources of systematic errors in the lattice evaluation and the estimate from 
uniisii, the agreement between the two is encouraging, especially for the isovector combination, where 
sea quark contributions should be small. 

The angular momentum carried by sea quarks is hardly known. It cannot be inferred form electro¬ 
magnetic form factors, and in lattice QCD it requires calculation of so-called “disconnected” graphs, 
which are affected by large statistical errors. Information on sea quarks can however be obtained 
from hard exclusive processes, especially from DVGS, where the control of theory is highest and a 
large number of spin asymmetries can be evaluated in the leading-twist approximation. HERMES 
has presented preliminary results for the spin asymmetry in DVGS on a transversely polarized proton 
|33) . Their comparison with a calculation m using J*? as input parameter in the ansatz for E‘> 
shows that this asymmetry is indeed sensitive to in the kinematics of the experiment (recall that 
in Compton scattering on the proton u and d quark distributions have relative weight 4 : 1 due to 
the squared quark charges). A measurement of DVGS on the neutron at Jefferson Lab Hall A |35j is 
currently being analyzed and should provide information about A number of ongoing measure¬ 
ments at Jefferson Lab and DESY can provide information on the unpolarized distributions H'^ and 

. Future experimental prospects for this field are a proposed run of COMPASS at CERN with a 
detector for recoiling protons ESI, the planned upgrade of Jefferson Lab to 12 GeV, and eventually 
an electron-proton collider eRHIC/EIC. 

5 Conclusions 

In this talk I have reviewed progress in selected areas connected with generalized parton distributions. 
There are many interesting aspects discussed in the recent literature which I could not cover for reasons 
of time. Among these are studies of GPDs in dynamical models, GPDs of nuclei, production of exotic 
mesons, GPDs for hadron-photon and baryon-meson transitions, hard exclusive processes at large 
s and t, and possibilities to study GPDs in neutrino-nucleon scattering with future high-intensity 
neutrino beams. 

To summarize, there has been important technical progress in the description of hard exclusive 
processes, with full NLO results in q;<j available for most relevant channels, partial NNLO results for 
Gompton scattering, and a better understanding of the scale evolution of GPDs. These advances 
remain to be fully implemented in the analysis of data, but existing studies have in particular shown 
the high sensitivity of vector meson production to the generalized gluon distribution, even at the 
moderate x accessible in fixed-target experiments. 

Through the impact parameter representation, GPDs provide information on the spatial distri¬ 
bution of partons in a hadron, and a number of studies have turned this concept into quantitative 
information. Phenomenological analysis of elastic form factors, as well as results from lattice QCD 
show that the average impact parameter of valence quarks strongly decreases with their momentum 
fraction in the proton. The proton helicity-flip distribution has connections with two crucial as¬ 
pects of spin physics: transverse polarization effects and the orbital angular momentum carried 
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by quarks in the nucleon. First steps have been taken towards a quantitative understanding of 
The present picture—suggested both by lattice results and model-dependent analysis of experimental 
data—is that for individual quark flavors may be substantial, but that when summed over flavors 
the orbital angular momentum carried by valence quarks contributes little to the nucleon spin. The 
proton spin puzzle thus remains a puzzle. 
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